- Author:
- Glucose Uptake <nafs080@aucklanduni.ac.nz>
- Date:
- 2018-03-29 11:43:57+13:00
- Desc:
- Modular_Version of glucose uptake model in CellML
- Permanent Source URI:
- http://models.cellml.org/workspace/572/rawfile/005feb6e7aaf2994df6bb8eadc6e07a310866f77/Composite Model(01.02.2018)(No Condition).cellml
<?xml version='1.0'?>
<model name="Nima_2018" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<!-- electroneutrality maintained also apical and basolateral concentration considered as constant without any equations
Basolateral concentration are constant and they are same as high glucose DMEM cell culture
No conditon for glucose concentration in the lumen
Output of this model then fitted to experimental data-->
<units name="per_second">
<unit exponent="-1" units="second"/>
</units>
<units name="per_second2">
<unit exponent="-1" units="second"/>
</units>
<!-- def unit per_second5 as
unit second {expo: -5};
enddef;-->
<units name="volt_per_second">
<unit units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_meter">
<unit exponent="-1" units="meter"/>
</units>
<units name="per_volt_per_second">
<unit exponent="-1" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt4_per_second">
<unit exponent="-4" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt3_per_second">
<unit exponent="-3" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt2_per_second">
<unit exponent="-2" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_mol">
<unit exponent="-1" units="mole"/>
</units>
<units name="umol">
<unit prefix="micro" units="mole"/>
</units>
<units name="per_umol">
<unit exponent="-1" units="umol"/>
</units>
<units name="M">
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="per_M">
<unit exponent="-1" units="M"/>
</units>
<units name="m2">
<unit exponent="2" units="metre"/>
</units>
<units name="cm">
<unit prefix="centi" units="metre"/>
</units>
<units name="um2">
<unit exponent="2" prefix="micro" units="metre"/>
</units>
<units name="um">
<unit prefix="micro" units="meter"/>
</units>
<units name="um3">
<unit exponent="3" prefix="micro" units="metre"/>
</units>
<units name="m3">
<unit exponent="3" units="metre"/>
</units>
<!-- def unit litre as
unit m3 {pref: milli};
enddef;-->
<units name="mol_per_um2">
<unit units="mole"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="mM">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="mV">
<unit prefix="milli" units="volt"/>
</units>
<units name="J_per_K_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="joule"/>
<unit exponent="-1" units="kelvin"/>
</units>
<units name="pA_per_cm2">
<unit prefix="pico" units="ampere"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="uA">
<unit prefix="micro" units="ampere"/>
</units>
<units name="uA_per_um2">
<unit units="uA"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="uA_per_A">
<unit units="uA"/>
<unit exponent="-1" units="ampere"/>
</units>
<units name="umol_per_mol">
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="C_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="coulomb"/>
</units>
<units name="per_M_per_second">
<unit exponent="-1" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M2_per_second">
<unit exponent="-2" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M4_per_second2">
<unit exponent="-4" units="M"/>
<unit exponent="-2" units="second"/>
</units>
<units name="per_M4_per_second">
<unit exponent="-4" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uM">
<unit prefix="micro" units="M"/>
</units>
<units name="M2">
<unit exponent="2" units="M"/>
</units>
<units name="M3">
<unit exponent="3" units="M"/>
</units>
<units name="M3_per_second">
<unit exponent="3" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M3_per_second5">
<unit exponent="-3" units="M"/>
<unit exponent="-5" units="second"/>
</units>
<units name="M_per_second">
<unit units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="cm2">
<unit exponent="2" prefix="centi" units="metre"/>
</units>
<units name="cm2_per_s">
<unit units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="mole_per_m2">
<unit units="mole"/>
<unit exponent="-1" units="m2"/>
</units>
<!-- def unit uM_per_cm2_per_s as
unit uM;
unit cm2 {expo: -1};
unit second {expo: -1};
enddef;-->
<units name="Farad">
<unit exponent="-1" units="volt"/>
<unit units="coulomb"/>
</units>
<units name="uF">
<unit prefix="micro" units="Farad"/>
</units>
<units name="uF_per_cm2">
<unit prefix="micro" units="Farad"/>
<unit exponent="-1" units="cm2"/>
</units>
<units name="umol_per_cm2_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="umol_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uA_per_cm2">
<unit units="uA"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="m_per_s">
<unit units="meter"/>
<unit exponent="-1" units="second"/>
</units>
<units name="m3_mole_per_litre_per_umol">
<unit units="m3"/>
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
<unit exponent="-1" units="umol"/>
</units>
<units name="litre_umol_per_mol_per_m3">
<unit units="litre"/>
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
<unit exponent="-1" units="m3"/>
</units>
<units name="S">
<unit units="siemens"/>
</units>
<units name="S_per_m2">
<unit units="siemens"/>
<unit exponent="-1" units="m2"/>
</units>
<units name="per_m2">
<unit exponent="-1" units="m2"/>
</units>
<units name="cm_per_s">
<unit units="cm"/>
<unit exponent="-1" units="second"/>
</units>
<component name="parameters">
<variable initial_value="-1" name="z_Cl" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_Na" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_k" public_interface="out" units="dimensionless"/>
<variable initial_value="96485" name="F" public_interface="out" units="C_per_mol"/>
<variable initial_value="8.314" name="R" public_interface="out" units="J_per_K_per_mol"/>
<variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
<variable initial_value="1e-5" name="capacitance" public_interface="out" units="uF"/>
<variable initial_value="8.4e-10" name="A_Apical" public_interface="out" units="m2"/>
<variable initial_value="5e-10" name="A_Basol" public_interface="out" units="m2"/>
<variable initial_value="5e-16" name="v_cell" public_interface="out" units="m3"/>
<variable initial_value="0" name="P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="4.5e-8" name="P_P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="1e-6" name="P_P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="3e-8" name="P_P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Cl" public_interface="out" units="m_per_s"/>
<variable name="time" public_interface="out" units="second"/>
<variable initial_value="14000" name="k0_12" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="300" name="k0_21" public_interface="out" units="per_second"/>
<variable initial_value="15000" name="k0_23" public_interface="out" units="per_M_per_second"/>
<variable initial_value="20" name="k0_32" public_interface="out" units="per_second"/>
<variable initial_value="0.1" name="k0_25" public_interface="out" units="per_second"/>
<variable initial_value="0.01" name="k0_52" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_34" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_43" public_interface="out" units="per_second"/>
<variable initial_value="6" name="k0_45" public_interface="out" units="per_second"/>
<variable initial_value="450" name="k0_54" public_interface="out" units="per_M_per_second"/>
<variable initial_value="10" name="k0_56" public_interface="out" units="per_second"/>
<variable initial_value="2000" name="k0_65" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="25" name="k0_61" public_interface="out" units="per_second"/>
<variable initial_value="600" name="k0_16" public_interface="out" units="per_second"/>
<variable initial_value="0.7" name="delta" public_interface="out" units="dimensionless"/>
<variable initial_value="0.3" name="alpha_p" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="alpha_pp" public_interface="out" units="dimensionless"/>
<!--ar N_C: dimensionless {init: 6e10};
var N_Avo: per_mol {init: 6.022e23};-->
<!--var C_T: umol {pub: out};-->
<variable initial_value="2" name="n" public_interface="out" units="dimensionless"/>
<variable initial_value="-2" name="z_c" public_interface="out" units="dimensionless"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="mu" public_interface="out" units="dimensionless"/>
<!--alpha_p = 1{dimensionless}-delta-alpha_pp;-->
<!--C_T = 1e6{umol_per_mol}*N_C/N_Avo;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Apical_concentrations">
<variable initial_value="105e-3" name="Na_m" public_interface="out" units="M"/>
<variable initial_value="25e-3" name="K_m" public_interface="out" units="M"/>
<variable initial_value="50e-3" name="glucose_m" public_interface="out" units="M"/>
<variable initial_value="106e-3" name="Cl_m" public_interface="out" units="M"/>
<!-- var k: per_second {init: 0.2, pub: out};
var glucose_0: M {init: 20e-3, pub: out};
var glucose_f: M {init: 0.1e-3, pub: out};
var time: second {pub: in};-->
<!-- glucose_m = sel
case (time >= 0{second}) and (time < 60{second}):
0.1e-3{M};
case time >= 60{second}:
(glucose_0-glucose_f)*exp(-k*(time-60{second}))+glucose_f;
endsel;-->
</component>
<component name="Apical_voltage">
<variable initial_value="-46e-3" name="v_mc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<variable name="I_D_Na" public_interface="in" units="uA"/>
<variable name="I_D_K" public_interface="in" units="uA"/>
<variable name="I_D_Cl" public_interface="in" units="uA"/>
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NHE3_H" public_interface="in" units="uA"/>
<variable name="I_NHE3_Na" public_interface="in" units="uA"/>
<variable name="I_BK" public_interface="in" units="uA"/>
<!-- var I_AE1_HCO3: uA {pub: in};
var I_AE1_Cl: uA {pub: in};-->
<variable name="I_ENaC" public_interface="in" units="uA"/>
<!--var I_NaCl: uA {pub: in};-->
<variable name="I_Na_SGLT" public_interface="in" units="uA"/>
<variable name="I_CFTR" public_interface="in" units="uA"/>
<!--var I_CFTR_HCO3: uA {pub: in};-->
<!--var I_Gl_SGLT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_mc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_D_Na</ci>
<ci>I_D_K</ci>
</apply>
<ci>I_D_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_Na_SGLT</ci>
</apply>
<ci>I_BK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_CFTR</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_ENaC</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_D_P_Na</ci>
<ci>I_D_P_K</ci>
</apply>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="Cell_concentration">
<variable initial_value="47e-3" name="Na_i" public_interface="out" units="M"/>
<!--var Na_m: M {pub: in};
var Na_s: M {pub: in};-->
<variable initial_value="11e-3" name="glucose_i" public_interface="out" units="M"/>
<variable initial_value="129e-3" name="K_i" public_interface="out" units="M"/>
<variable initial_value="60e-6" name="Ca_i" public_interface="out" units="M"/>
<variable initial_value="57e-3" name="Cl_i" public_interface="out" units="M"/>
<variable name="HCO3_i" public_interface="out" units="M"/>
<variable name="H_int" public_interface="out" units="M"/>
<variable name="HCO3_m" public_interface="out" units="M"/>
<variable name="HCO3_s" public_interface="out" units="M"/>
<variable name="H_Ext" public_interface="out" units="M"/>
<variable name="H_s" public_interface="out" units="M"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="J_Na_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Gl_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="in" units="umol_per_s"/>
<!--var J_NaCl: umol_per_s {pub: in};-->
<variable name="J_NHE3_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_A_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_D_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_AE1_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_AE1_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_CFTR" public_interface="in" units="umol_per_s"/>
<variable name="J_CFTR_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_ENaC" public_interface="in" units="umol_per_s"/>
<variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
<variable name="J_IK" public_interface="in" units="umol_per_s"/>
<variable name="J_BK" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
<!-- var J_buf_C: M_per_second {pub: in};
var J_buf_A: M_per_second {pub: in};-->
<variable name="time" public_interface="in" units="second"/>
<variable initial_value="7.4" name="pH_int" public_interface="out" units="per_second"/>
<variable initial_value="7.4" name="pH_Ext" public_interface="out" units="per_second"/>
<variable initial_value="7.4" name="pH_s" public_interface="out" units="per_second"/>
<variable initial_value="1.68e6" name="rho" public_interface="out" units="per_meter"/>
<variable name="beta" public_interface="out" units="M"/>
<variable initial_value="6.1" name="p_K" public_interface="out" units="dimensionless"/>
<variable initial_value="40" name="P_CO2" public_interface="out" units="dimensionless"/>
<variable initial_value="0.03e-3" name="s" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="theta_1" units="dimensionless"/>
<variable initial_value="1" name="theta_2" units="dimensionless"/>
<variable initial_value="1" name="theta_3" units="dimensionless"/>
<variable initial_value="1" name="theta_4" units="dimensionless"/>
<variable initial_value="1" name="theta_5" units="dimensionless"/>
<variable initial_value="1" name="theta_6" units="dimensionless"/>
<variable initial_value="1" name="theta_7" units="dimensionless"/>
<variable initial_value="1" name="theta_8" units="dimensionless"/>
<variable initial_value="1" name="theta_9" units="dimensionless"/>
<variable initial_value="1" name="theta_10" units="dimensionless"/>
<variable initial_value="1" name="theta_11" units="dimensionless"/>
<variable initial_value="1" name="theta_12" units="dimensionless"/>
<variable initial_value="1" name="theta_13" units="dimensionless"/>
<variable initial_value="1" name="theta_14" units="dimensionless"/>
<variable initial_value="1" name="theta_15" units="dimensionless"/>
<variable initial_value="1" name="theta_16" units="dimensionless"/>
<variable initial_value="1" name="theta_17" units="dimensionless"/>
<variable initial_value="1" name="theta_18" units="dimensionless"/>
<variable initial_value="1" name="theta_19" units="dimensionless"/>
<variable initial_value="1" name="theta_20" units="dimensionless"/>
<variable initial_value="1" name="theta_21" units="dimensionless"/>
<variable initial_value="1" name="theta_22" units="dimensionless"/>
<variable initial_value="1" name="theta_23" units="dimensionless"/>
<variable initial_value="1" name="theta_24" units="dimensionless"/>
<variable initial_value="1" name="theta_25" units="dimensionless"/>
<variable initial_value="1" name="theta_26" units="dimensionless"/>
<variable initial_value="1" name="theta_27" units="dimensionless"/>
<variable name="K_Concentration" public_interface="out" units="M"/>
<variable name="Na_Concentration" public_interface="out" units="M"/>
<variable name="Cl_Concentration" public_interface="out" units="M"/>
<!--var Concentration_Na: M {pub: out};-->
<!-- Concentration_Na = Na_m+Na_i+Na_s;-->
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Na_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_1</ci>
<ci>J_Na_SGLT</ci>
</apply>
<apply>
<times/>
<ci>theta_2</ci>
<ci>J_Na_NaK</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_15</ci>
<ci>J_ENaC</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_16</ci>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_19</ci>
<ci>J_NBC_Na</ci>
</apply>
</apply>
<ci>J_D_P_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_K</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_Bl_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0</cn>
<ci>J_IK</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Cl_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_11</ci>
<ci>J_CFTR</ci>
</apply>
<apply>
<times/>
<ci>theta_18</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_12</ci>
<ci>J_D_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_13</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2.3</cn>
<ci>HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_1</ci>
<ci>J_Na_SGLT</ci>
</apply>
<apply>
<times/>
<ci>theta_2</ci>
<ci>J_Na_NaK</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_15</ci>
<ci>J_ENaC</ci>
</apply>
<apply>
<times/>
<ci>theta_16</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_19</ci>
<ci>J_NBC_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_4</ci>
<ci>J_D_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_5</ci>
<ci>J_D_Bl_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>glucose_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_26</ci>
<ci>J_A_GLUT</ci>
</apply>
<apply>
<times/>
<ci>theta_6</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_7</ci>
<ci>J_GLUT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_8</ci>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<ci>theta_17</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_9</ci>
<ci>J_D_K</ci>
</apply>
<apply>
<times/>
<ci>theta_10</ci>
<ci>J_D_Bl_K</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_27</ci>
<ci>J_IK</ci>
</apply>
</apply>
<ci>J_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cl_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_11</ci>
<apply>
<minus/>
<ci>J_CFTR</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_18</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_12</ci>
<ci>J_D_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_13</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<!--ode(HCO3_i, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*theta_20*J_NBC_HCO3+theta_25*J_AE1_HCO3+theta_21*J_buf_C;-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pH_int</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>rho</ci>
</apply>
<apply>
<times/>
<ci>beta</ci>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_20</ci>
<ci>J_NBC_HCO3</ci>
</apply>
<apply>
<times/>
<ci>theta_22</ci>
<ci>J_NHE3_H</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_25</ci>
<ci>J_AE1_HCO3</ci>
</apply>
</apply>
</apply>
</apply>
<!--ode(pH_Ext, time) = 1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*(theta_22*J_NHE3_H-theta_25*J_AE1_HCO3-theta_26*J_CFTR_HCO3);-->
<!--ode(pH_s, time) = -1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*theta_20*J_NBC_HCO3;-->
<!--ode(HCO3_m, time) = theta_24*-J_AE1_HCO3+theta_23*J_buf_A;-->
<apply>
<eq/>
<ci>H_int</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_i</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_Ext</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_Ext</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_s</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_m</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_Ext</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_s</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_s</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Basolateral_concentrations">
<variable initial_value="155e-3" name="Na_s" public_interface="out" units="M"/>
<!--var glucose_s: M {init: 10e-3, pub: out};-->
<variable initial_value="119e-3" name="Cl_s" public_interface="out" units="M"/>
<variable initial_value="5.4e-3" name="K_s" public_interface="out" units="M"/>
<variable initial_value="25e-3" name="glucose_s" public_interface="out" units="M"/>
<!--var time: second {pub: out};-->
</component>
<component name="Basol_voltage">
<variable initial_value="-50e-3" name="v_sc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<variable name="I_D_Bl_Na" public_interface="in" units="uA"/>
<variable name="I_D_Bl_K" public_interface="in" units="uA"/>
<variable name="I_D_Bl_Cl" public_interface="in" units="uA"/>
<variable name="I_Na_NaK" public_interface="in" units="uA"/>
<variable name="I_K_NaK" public_interface="in" units="uA"/>
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NBC_Na" public_interface="in" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="in" units="uA"/>
<variable name="I_IK" public_interface="in" units="uA"/>
<!--var I_GLUT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_sc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>I_D_Bl_Na</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_Bl_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_Bl_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<minus/>
<ci>I_Na_NaK</ci>
</apply>
<ci>I_K_NaK</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_HCO3</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_IK</ci>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="rate_constants">
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable name="k_12" public_interface="out" units="per_second"/>
<variable name="ks_12" public_interface="out" units="per_M2_per_second"/>
<variable name="k_21" public_interface="out" units="per_second"/>
<!--var k_23: per_second {pub: out};-->
<variable name="k_32" public_interface="out" units="per_second"/>
<variable name="k_25" public_interface="out" units="per_second"/>
<variable name="k_52" public_interface="out" units="per_second"/>
<variable name="k_34" public_interface="out" units="per_second"/>
<variable name="k_43" public_interface="out" units="per_second"/>
<variable name="k_45" public_interface="out" units="per_second"/>
<variable name="k_54" public_interface="out" units="per_second"/>
<variable name="k_56" public_interface="out" units="per_second"/>
<variable name="k_65" public_interface="out" units="per_second"/>
<variable name="k_61" public_interface="out" units="per_second"/>
<variable name="k_16" public_interface="out" units="per_second"/>
<variable name="k0_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_21" public_interface="in" units="per_second"/>
<!--var k0_23: per_M_per_second {pub: in};-->
<variable name="k0_32" public_interface="in" units="per_second"/>
<variable name="k0_25" public_interface="in" units="per_second"/>
<variable name="k0_52" public_interface="in" units="per_second"/>
<variable name="k0_34" public_interface="in" units="per_second"/>
<variable name="k0_43" public_interface="in" units="per_second"/>
<variable name="k0_45" public_interface="in" units="per_second"/>
<variable name="k0_54" public_interface="in" units="per_M_per_second"/>
<!--var k0_25: per_second {pub: in};
var k0_52: per_second {pub: in};-->
<variable name="k0_56" public_interface="in" units="per_second"/>
<variable name="k0_65" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_61" public_interface="in" units="per_second"/>
<variable name="k0_16" public_interface="in" units="per_second"/>
<!--var k0_54_temp: per_M_per_second {pub: out};-->
<!--var k_52_temp: per_second {pub: out};-->
<variable name="delta" public_interface="in" units="dimensionless"/>
<variable name="n" public_interface="in" units="dimensionless"/>
<variable name="z_c" public_interface="in" units="dimensionless"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="alpha_p" public_interface="in" units="dimensionless"/>
<!--var alpha_pp: dimensionless {pub: in};-->
<!--var C_T: umol {pub: in};-->
<variable name="mu" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ks_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>n</ci>
</apply>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<!-- k_12 = ks_12*pow(Na_m, n);-->
<apply>
<eq/>
<ci>k_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha_p</ci>
</apply>
<ci>mu</ci>
</apply>
</apply>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_21</ci>
<apply>
<times/>
<ci>k0_21</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>n</ci>
<ci>z_Na</ci>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<!--k_23 = k0_23*glucose_m;-->
<apply>
<eq/>
<ci>k_32</ci>
<ci>k0_32</ci>
</apply>
<apply>
<eq/>
<ci>k_34</ci>
<ci>k0_34</ci>
</apply>
<apply>
<eq/>
<ci>k_43</ci>
<ci>k0_43</ci>
</apply>
<apply>
<eq/>
<ci>k_45</ci>
<ci>k0_45</ci>
</apply>
<apply>
<eq/>
<ci>k_54</ci>
<apply>
<times/>
<ci>k0_54</ci>
<ci>glucose_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_25</ci>
<ci>k0_25</ci>
</apply>
<!--k_52_temp = k0_52*exp((z_c+n)*delta*mu/2{dimensionless});-->
<apply>
<eq/>
<ci>k_56</ci>
<ci>k0_56</ci>
</apply>
<apply>
<eq/>
<ci>k_65</ci>
<apply>
<times/>
<ci>k0_65</ci>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_61</ci>
<apply>
<times/>
<ci>k0_61</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_c</ci>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_16</ci>
<apply>
<times/>
<ci>k0_16</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>z_c</ci>
</apply>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_52</ci>
<ci>k0_52</ci>
</apply>
</math>
<!--k0_54_temp = k0_23*k_34*k_45/(k_43*k_32);-->
</component>
<component name="phenomonological_constants">
<variable name="epsilon" public_interface="out" units="per_second"/>
<variable name="lambda" public_interface="out" units="per_M3_per_second5"/>
<variable name="chi" public_interface="out" units="M"/>
<variable name="alpha" public_interface="out" units="M3"/>
<variable name="beta" public_interface="out" units="M2"/>
<variable name="gamma" public_interface="out" units="M3_per_second"/>
<variable name="phi" public_interface="out" units="M_per_second"/>
<variable name="k_12" public_interface="in" units="per_second"/>
<variable name="k_21" public_interface="in" units="per_second"/>
<!--var k_23: per_second {pub: in};-->
<variable name="k_32" public_interface="in" units="per_second"/>
<variable name="k_34" public_interface="in" units="per_second"/>
<variable name="k_43" public_interface="in" units="per_second"/>
<variable name="k_45" public_interface="in" units="per_second"/>
<variable name="k_54" public_interface="in" units="per_second"/>
<variable name="k_25" public_interface="in" units="per_second"/>
<variable name="k_52" public_interface="in" units="per_second"/>
<variable name="k_56" public_interface="in" units="per_second"/>
<variable name="k_65" public_interface="in" units="per_second"/>
<variable name="k_61" public_interface="in" units="per_second"/>
<variable name="k_16" public_interface="in" units="per_second"/>
<variable name="ks_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_23" public_interface="in" units="per_M_per_second"/>
<!--var C_T: umol {pub: in};-->
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable initial_value="22e7" name="n_SGLT" public_interface="out" units="dimensionless"/>
<variable name="J_Na_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="I_Na_SGLT" public_interface="out" units="uA"/>
<variable name="J_Gl_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="J_GL_SGLT1" public_interface="out" units="umol_per_s"/>
<!--var I_Gl_SGLT: uA {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>lambda</ci>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>chi</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>k0_23</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_65</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>phi</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k_56</ci>
<ci>k_61</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_SGLT</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>n_SGLT</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>epsilon</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>phi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>gamma</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>alpha</ci>
<apply>
<times/>
<ci>beta</ci>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>chi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_SGLT</ci>
<apply>
<times/>
<ci>J_Na_SGLT</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Gl_SGLT</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>J_Na_SGLT</ci>
</apply>
</apply>
<!--to change the unit to picomol/hour-->
<apply>
<eq/>
<ci>J_GL_SGLT1</ci>
<apply>
<times/>
<ci>J_Gl_SGLT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
</math>
<!--I_Gl_SGLT = -J_Gl_SGLT*F;-->
</component>
<component name="CFTR">
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="V_A_Cl" public_interface="out" units="volt"/>
<variable name="V_A_HCO3" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_CFTR" public_interface="out" units="umol_per_s"/>
<variable name="J_CFTR1" public_interface="out" units="umol_per_s"/>
<variable name="J_CFTR_HCO3" public_interface="out" units="umol_per_s"/>
<variable name="I_CFTR" public_interface="out" units="uA"/>
<variable name="I_CFTR_HCO3" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="10" name="G_CFTR" public_interface="out" units="S"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_m" public_interface="in" units="M"/>
<variable initial_value="-1" name="z_HCO3" public_interface="out" units="dimensionless"/>
<!-- alpha is multipled by area in order to have a dimensionless unit-->
<variable initial_value="0.1" name="alpha_CFTR" public_interface="out" units="per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_CFTR</ci>
<apply>
<divide/>
<ci>I_CFTR</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_CFTR1</ci>
<apply>
<times/>
<ci>J_CFTR</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CFTR</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_CFTR</ci>
<ci>alpha_CFTR</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Cl</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_HCO3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>HCO3_m</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CFTR_HCO3</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">0.25<sep/>6</cn>
<ci>G_CFTR</ci>
<ci>alpha_CFTR</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_HCO3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_CFTR_HCO3</ci>
<apply>
<divide/>
<ci>I_CFTR_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="ENaC">
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="V_A_Na" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_ENaC" public_interface="out" units="umol_per_s"/>
<variable name="J_ENaC1" public_interface="out" units="umol_per_s"/>
<variable name="I_ENaC" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="11" name="G_ENaC" public_interface="out" units="S"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<variable initial_value="0.1" name="alpha_ENaC" public_interface="out" units="per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_ENaC</ci>
<apply>
<divide/>
<ci>I_ENaC</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_ENaC1</ci>
<apply>
<times/>
<ci>J_ENaC</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_ENaC</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_ENaC</ci>
<ci>alpha_ENaC</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaKCC">
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="S_C" public_interface="out" units="M"/>
<variable name="J_NaKCC" public_interface="out" units="umol_per_s"/>
<variable name="I_NaKCC" public_interface="out" units="uA"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="u_NaKCC" public_interface="out" units="per_second"/>
<variable initial_value="157.55" name="gamma_1" public_interface="out" units="per_second2"/>
<variable initial_value="2e7" name="gamma_2" public_interface="out" units="per_M4_per_second2"/>
<variable initial_value="1.03" name="gamma_3" public_interface="out" units="per_second"/>
<variable initial_value="1.38e6" name="gamma_4" public_interface="out" units="per_M4_per_second"/>
<variable initial_value="8e-9" name="alpha_NaKCC" public_interface="out" units="mole_per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>S_C</ci>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Cl_i</ci>
<ci>Cl_i</ci>
<ci>K_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_NaKCC</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>gamma_1</ci>
<apply>
<times/>
<ci>gamma_2</ci>
<ci>S_C</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>gamma_3</ci>
<apply>
<times/>
<ci>gamma_4</ci>
<ci>S_C</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaKCC</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>alpha_NaKCC</ci>
<ci>u_NaKCC</ci>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NaKCC</ci>
<apply>
<times/>
<ci>J_NaKCC</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="A_GLUT2">
<variable initial_value="1" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="11130" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="1" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="11" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_A_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="mole"/>
<variable initial_value="9e7" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_A_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<!-- D = sel
case (glucose_m > 0.02{M}):
(1{dimensionless}+glucose_i/K_S)*(K_32*glucose_m/K_S+K_41)+(1{dimensionless}+glucose_m/K_S)*(K_23*glucose_i/K_S+K_14);
case glucose_m < 0.02{M}:
1e100{per_second};
endsel;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_mole" type="e-notation">6.022<sep/>23</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_A_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-9</cn>
<ci>C_tot</ci>
</apply>
<apply>
<times/>
<ci>D</ci>
<ci>v_cell</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_A_GLUT1</ci>
<apply>
<times/>
<ci>J_A_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_m</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="GLUT2">
<variable initial_value="1" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="2000" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="0.726" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="12.1" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_s" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="mole"/>
<variable initial_value="8e7" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_mole" type="e-notation">6.022<sep/>23</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-9</cn>
<ci>C_tot</ci>
</apply>
<apply>
<times/>
<ci>D</ci>
<ci>v_cell</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_GLUT1</ci>
<apply>
<times/>
<ci>J_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_s</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaK_channel">
<variable name="time" public_interface="in" units="second"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="TO_NaK" public_interface="out" units="per_second"/>
<variable initial_value="5e6" name="n_NaK" public_interface="out" units="dimensionless"/>
<variable initial_value="10e-3" name="K_Na_i_half" public_interface="out" units="M"/>
<variable initial_value="1.5e-3" name="K_K_s_half" public_interface="out" units="M"/>
<variable name="I_Na_NaK" public_interface="out" units="uA"/>
<variable name="I_K_NaK" public_interface="out" units="uA"/>
<variable name="J_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="out" units="umol_per_s"/>
<!--var J_K_NaK: umol_per_s {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_K_NaK</ci>
<apply>
<times/>
<ci>J_K_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_NaK</ci>
<apply>
<times/>
<ci>J_Na_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>TO_NaK</ci>
<ci>n_NaK</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<cn cellml:units="dimensionless">1.13</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_Na_i_half</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1.3</cn>
<ci>K_s</ci>
</apply>
<apply>
<plus/>
<ci>K_K_s_half</ci>
<ci>K_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>TO_NaK</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_volt4_per_second" type="e-notation">5.46<sep/>+4</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt3_per_second" type="e-notation">2.43<sep/>+2</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt2_per_second" type="e-notation">2<sep/>+3</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt_per_second">160</cn>
<ci>v_sc</ci>
</apply>
<cn cellml:units="per_second">51</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>J_NaK</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_K_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaK</ci>
</apply>
</apply>
</math>
<!--J_Na_NaK = -108e8{dimensionless}*J_NaK;-->
<!--J_K_NaK = -72e8{dimensionless}*J_NaK;-->
</component>
<component name="NBC">
<variable name="time" public_interface="in" units="second"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable initial_value="10e-3" name="HCO3_s" public_interface="out" units="M"/>
<variable initial_value="0.01" name="g_NBC" public_interface="out" units="S_per_m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<!--kelvin {init: 300};-->
<variable initial_value="2" name="n" units="dimensionless"/>
<!--var eps: M {init: 1e-12};-->
<variable name="E_NBC" public_interface="out" units="volt"/>
<variable name="J_NBC_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="out" units="umol_per_s"/>
<variable name="I_NBC_Na" public_interface="out" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="out" units="uA"/>
<variable initial_value="1e-18" name="eps" units="M"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_NBC</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>n</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<times/>
<apply>
<divide/>
<ci>Na_s</ci>
<ci>Na_i</ci>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>HCO3_s</ci>
<ci>HCO3_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<ci>g_NBC</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>E_NBC</ci>
<ci>v_sc</ci>
</apply>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_HCO3</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3.0</cn>
<ci>J_NBC_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_Na</ci>
<apply>
<times/>
<ci>J_NBC_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_HCO3</ci>
<apply>
<times/>
<ci>J_NBC_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="NHE3">
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="H_Ext" public_interface="in" units="M"/>
<variable name="H_int" public_interface="in" units="M"/>
<variable initial_value="0.0" name="NH4_ext" public_interface="out" units="M"/>
<variable initial_value="0.0" name="NH4_int" public_interface="out" units="M"/>
<variable name="J_NHE3_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H1" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_NH4" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_Na_Max" public_interface="out" units="umol_per_s"/>
<variable name="I_NHE3_H" public_interface="out" units="uA"/>
<variable name="I_NHE3_Na" public_interface="out" units="uA"/>
<variable initial_value="1e-3" name="x_T" units="M"/>
<variable name="sigma" units="cm_per_s"/>
<variable name="P_Na" units="cm_per_s"/>
<variable name="P_H" units="cm_per_s"/>
<variable name="P_NH4" units="cm_per_s"/>
<variable initial_value="1.6e-3" name="P0_Na" units="cm_per_s"/>
<variable initial_value="0.48e-3" name="P0_H" units="cm_per_s"/>
<variable initial_value="1.6e-3" name="P0_NH4" units="cm_per_s"/>
<variable initial_value="30e-3" name="K_Na" units="M"/>
<variable initial_value="72e-9" name="K_H" units="M"/>
<variable initial_value="27e-3" name="K_NH4" units="M"/>
<variable initial_value="1.0e-9" name="K_I" units="M"/>
<variable initial_value="0.0" name="f_m" units="dimensionless"/>
<variable initial_value="2.0" name="f_M" units="dimensionless"/>
<variable name="alpha_ext_Na" units="dimensionless"/>
<variable name="alpha_int_Na" units="dimensionless"/>
<variable name="beta_ext_H" units="dimensionless"/>
<variable name="beta_int_H" units="dimensionless"/>
<variable name="gamma_ext_NH4" units="dimensionless"/>
<variable name="gamma_int_NH4" units="dimensionless"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<!-- var pH_Ext: dimensionless;
var pH_int: dimensionless;
var p_K: dimensionless {init: 6.1, pub: out};
var P_CO2: dimensionless {init: 40, pub: out};
var s: dimensionless {init: 0.03e-3, pub: out};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};-->
<!-- pH_Ext = p_K+log(HCO3_m/(s*P_CO2));
pH_int = p_K+log(HCO3_i/(s*P_CO2));-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_ext_Na</ci>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_int_Na</ci>
<apply>
<divide/>
<ci>Na_i</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext_H</ci>
<apply>
<divide/>
<ci>H_Ext</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int_H</ci>
<apply>
<divide/>
<ci>H_int</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext_NH4</ci>
<apply>
<divide/>
<ci>NH4_ext</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int_NH4</ci>
<apply>
<divide/>
<ci>NH4_int</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_Na</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_H</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_H</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_NH4</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_NH4</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_ext_Na</ci>
<ci>beta_ext_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_int_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_int_Na</ci>
<ci>beta_int_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_ext_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>7</cn>
</apply>
<ci>A_Apical</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>7</cn>
<ci>A_Apical</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_int_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H1</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_NH4</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>gamma_ext_NH4</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na_Max</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
<apply>
<plus/>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_H</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_Na</ci>
<apply>
<times/>
<ci>J_NHE3_Na</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="AE1">
<!--var A_Apical: m2 {pub: in};-->
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_m" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="J_AE1_HCO3" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_HCO31: umol_per_s {pub: out};-->
<variable name="J_AE1_Cl" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_Cl1: umol_per_s {pub: out};-->
<variable name="I_AE1_HCO3" public_interface="out" units="uA"/>
<variable name="I_AE1_Cl" public_interface="out" units="uA"/>
<variable name="J_HCO3_influx" public_interface="out" units="umol_per_s"/>
<variable name="J_Cl_influx" public_interface="out" units="umol_per_s"/>
<variable initial_value="198e-3" name="K_HCO3_m" units="M"/>
<variable initial_value="198e-3" name="K_HCO3_i" units="M"/>
<variable initial_value="50e-3" name="K_Cl_s" units="M"/>
<variable initial_value="50e-3" name="K_Cl_i" units="M"/>
<variable initial_value="1.247" name="P_HCO3_m" units="cm_per_s"/>
<variable initial_value="0.135" name="P_HCO3_i" units="cm_per_s"/>
<variable initial_value="0.562" name="P_Cl_s" units="cm_per_s"/>
<variable initial_value="0.061" name="P_Cl_i" units="cm_per_s"/>
<variable name="beta_ext" units="dimensionless"/>
<variable name="beta_int" units="dimensionless"/>
<variable name="gamma_ext" units="dimensionless"/>
<variable name="gamma_int" units="dimensionless"/>
<variable name="sigma" units="cm_per_s"/>
<variable initial_value="4e-6" name="x_Tmax" units="M"/>
<variable initial_value="172e-3" name="K_I" units="M"/>
<variable name="x_T" units="M"/>
<variable name="x_ext" units="M"/>
<variable name="x_int" units="M"/>
<variable name="Jo_bm" units="umol_per_s"/>
<variable name="Ji_bm" units="umol_per_s"/>
<variable name="Js_bm" units="umol_per_s"/>
<variable name="Jo_cm" units="umol_per_s"/>
<variable name="Ji_cm" units="umol_per_s"/>
<variable name="Js_cm" units="umol_per_s"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>x_T</ci>
<apply>
<divide/>
<ci>x_Tmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext</ci>
<apply>
<divide/>
<ci>HCO3_m</ci>
<ci>K_HCO3_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int</ci>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext</ci>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>K_Cl_s</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int</ci>
<apply>
<divide/>
<ci>Cl_i</ci>
<ci>K_Cl_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_ext</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_int</ci>
<ci>gamma_int</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_AE1_HCO3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>7</cn>
<ci>A_Apical</ci>
<ci>x_Tmax</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<!--J_AE1_HCO31 = J_AE1_HCO3*36e8{dimensionless};
J_AE1_Cl1 = -J_AE1_HCO31;-->
<apply>
<eq/>
<ci>J_AE1_Cl</ci>
<ci>J_AE1_HCO3</ci>
</apply>
<apply>
<eq/>
<ci>I_AE1_HCO3</ci>
<apply>
<times/>
<ci>J_AE1_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_AE1_Cl</ci>
<apply>
<times/>
<ci>J_AE1_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_ext</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_int</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_HCO3_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Cl_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_i</ci>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_m</ci>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>HCO3_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_i</ci>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_s</ci>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>Cl_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<!-- def comp Hydration as
var time: second {pub: in};
var v_cell: m3 {pub: in};
var CO2_i: M {init: 1e-3, pub: out};
var CO2_s: M {init: 1.2e-3, pub: out};
var CO2_m: M {init: 1.2e-3, pub: out};
var H_int: M {pub: in};
var H_Ext: M {pub: in};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};
var p_CO2: per_second {init: 5e-6, pub: out};
var k_f: per_second {init: 0.03, pub: out};
var k_r: per_M_per_second {init: 20e-3, pub: out};
var J_CDF_A: umol_per_s {pub: out};
var J_CDF_B: umol_per_s {pub: out};
var J_buf_C: umol_per_s {pub: out};
var J_buf_C1: umol_per_s {pub: out};
var J_buf_A: umol_per_s {pub: out};
var J_D_P_CO2: umol_per_s {pub: in};-->
<!-- J_CDF_A = -1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_m);
J_CDF_B = 1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_s);
J_buf_C = 1e9{dimensionless}*v_cell*(k_f*CO2_i-k_r*HCO3_i*H_int);
J_buf_C1 = J_buf_C*36e8{dimensionless};
J_buf_A = 1e9{dimensionless}*v_cell*(k_f*CO2_m-k_r*HCO3_m*H_Ext);
ode(CO2_i, time) = -J_buf_C+J_CDF_A+J_CDF_B;
ode(CO2_m, time) = -J_buf_A-J_CDF_A-0{dimensionless}*J_D_P_CO2;
enddef;-->
<component name="IK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_IK" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="J_IK" public_interface="out" units="umol_per_s"/>
<variable name="I_IK" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="50" name="G_IK" public_interface="out" units="S"/>
<variable name="P_IK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.3e-3" name="K_s_IK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<variable initial_value="1" name="alpha_IK" public_interface="out" units="per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_IK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_s_IK</ci>
<ci>Ca_i</ci>
</apply>
<ci>eta</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_IK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_s</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_IK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
<ci>G_IK</ci>
<ci>alpha_IK</ci>
<ci>P_IK</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>V_A_IK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_IK</ci>
<apply>
<divide/>
<ci>I_IK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="BK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_BK" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_BK" public_interface="out" units="umol_per_s"/>
<variable name="I_BK" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="30" name="G_BK" public_interface="out" units="S"/>
<variable name="P_BK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.3e-3" name="K_s_BK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<variable initial_value="0.1" name="alpha_BK" public_interface="out" units="per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_BK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_s_BK</ci>
<ci>Ca_i</ci>
</apply>
<ci>eta</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_BK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_m</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_BK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
<ci>G_BK</ci>
<ci>alpha_BK</ci>
<ci>P_BK</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_BK</ci>
<apply>
<divide/>
<ci>I_BK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="Apical_Diffusion">
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_k" public_interface="in" units="m_per_s"/>
<variable name="P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_D_NA" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Na" public_interface="out" units="uA"/>
<variable name="J_D_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_K" public_interface="out" units="uA"/>
<variable name="J_D_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Apical" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Na_m</ci>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_NA</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
<ci>J_D_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Na</ci>
<apply>
<times/>
<ci>J_D_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>K_m</ci>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_K</ci>
<apply>
<times/>
<ci>J_D_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cl_m</ci>
<apply>
<times/>
<ci>Cl_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Cl</ci>
<apply>
<times/>
<ci>J_D_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Apical</ci>
<apply>
<plus/>
<ci>J_D_Na</ci>
<ci>J_D_K</ci>
<ci>J_D_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Basolateral_Diffusion">
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Basol" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Bl_Na" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_k" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Bl_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Na" public_interface="out" units="uA"/>
<variable name="J_D_Bl_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_K" public_interface="out" units="uA"/>
<variable name="J_D_Bl_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Basolateral" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Bl_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Na_s</ci>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Na</ci>
<apply>
<times/>
<ci>J_D_Bl_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>K_s</ci>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_K</ci>
<apply>
<times/>
<ci>J_D_Bl_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cl_s</ci>
<apply>
<times/>
<ci>Cl_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Cl</ci>
<apply>
<times/>
<ci>J_D_Bl_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Basolateral</ci>
<apply>
<plus/>
<ci>J_D_Bl_Na</ci>
<ci>J_D_Bl_K</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Paracellular_voltage">
<variable name="v_ms" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<!-- var time: second {pub: in};
var capacitance: uF {pub: in};
var I_D_P_Na: uA {pub: in};
var I_D_P_K: uA {pub: in};
var I_D_P_Cl: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ms</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>v_sc</ci>
</apply>
</apply>
</math>
<!--ode(v_ms, time) = (I_D_P_Na+I_D_P_K+I_D_P_Cl)/capacitance;-->
</component>
<component name="Paracellular_Diffusion">
<variable name="v_ms" public_interface="in" units="volt"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<!-- var CO2_m: M {pub: in};
var CO2_i: M {pub: in};-->
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_P_k" public_interface="in" units="m_per_s"/>
<variable name="P_P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_P_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Na" public_interface="out" units="uA"/>
<variable name="J_D_P_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_K" public_interface="out" units="uA"/>
<variable name="J_D_P_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Paracellular" public_interface="out" units="umol_per_s"/>
<!-- var I_D_P_CO2: uA {pub: out};
var J_D_P_CO2: umol_per_s {pub: out};
var G_P_CO2: S_per_m2 {init: 7, pub: out};
var v_P_CO2: volt {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_P_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Na_m</ci>
<apply>
<times/>
<ci>Na_s</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Na</ci>
<apply>
<times/>
<ci>J_D_P_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>K_m</ci>
<apply>
<times/>
<ci>K_s</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_K</ci>
<apply>
<times/>
<ci>J_D_P_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cl_m</ci>
<apply>
<times/>
<ci>Cl_s</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_ms</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Cl</ci>
<apply>
<times/>
<ci>J_D_P_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<!-- v_P_CO2 = R*T/F*ln(CO2_m/CO2_i);
I_D_P_CO2 = 1e6{dimensionless}*G_P_CO2*A_Apical*(v_ms-v_P_CO2);
J_D_P_CO2 = I_D_P_CO2/F;-->
<apply>
<eq/>
<ci>J_D_Paracellular</ci>
<apply>
<plus/>
<ci>J_D_P_Na</ci>
<ci>J_D_P_K</ci>
<ci>J_D_P_Cl</ci>
</apply>
</apply>
</math>
</component>
<!-- def comp interstitium_concentration as
var glucose_s: M {init: 25e-3, pub: out};-->
<!-- //var J_GLUT: umol_per_s {pub: in};
//var J_G_Blood: umol_per_s {pub: in};
//var v_cell: m3 {pub: in};
//var time: second {pub: in};-->
<!-- //ode(glucose_s, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*(-J_GLUT-J_G_Blood);
enddef;-->
<!-- def comp Blood_concentration as
var glucose_s: M {pub: in};
var glucose_B: M {init: 12e-3, pub: out};
var J_G_Blood: umol_per_s {pub: out};
var K: dimensionless {init: 0.000001, pub: out};-->
<!-- //var v_cell: m3 {pub: in};-->
<!-- var time: second {pub: in};-->
<!-- J_G_Blood = K*(glucose_s-glucose_B);
ode(glucose_B, time) = 0{dimensionless};
enddef;-->
<connection>
<map_components component_1="rate_constants" component_2="phenomonological_constants"/>
<map_variables variable_1="ks_12" variable_2="ks_12"/>
<map_variables variable_1="k_12" variable_2="k_12"/>
<map_variables variable_1="k_21" variable_2="k_21"/>
<!--vars k_23 and k_23;-->
<map_variables variable_1="k_32" variable_2="k_32"/>
<map_variables variable_1="k_25" variable_2="k_25"/>
<map_variables variable_1="k_52" variable_2="k_52"/>
<map_variables variable_1="k_34" variable_2="k_34"/>
<map_variables variable_1="k_43" variable_2="k_43"/>
<map_variables variable_1="k_45" variable_2="k_45"/>
<map_variables variable_1="k_54" variable_2="k_54"/>
<map_variables variable_1="k_56" variable_2="k_56"/>
<map_variables variable_1="k_65" variable_2="k_65"/>
<map_variables variable_1="k_61" variable_2="k_61"/>
<map_variables variable_1="k_16" variable_2="k_16"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Cell_concentration"/>
<map_variables variable_1="J_Na_SGLT" variable_2="J_Na_SGLT"/>
<map_variables variable_1="J_Gl_SGLT" variable_2="J_Gl_SGLT"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="parameters"/>
<map_variables variable_1="k0_23" variable_2="k0_23"/>
<!--vars C_T and C_T;-->
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_voltage"/>
<map_variables variable_1="I_Na_SGLT" variable_2="I_Na_SGLT"/>
<!--vars I_Gl_SGLT and I_Gl_SGLT;-->
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="parameters"/>
<map_variables variable_1="k0_12" variable_2="k0_12"/>
<map_variables variable_1="k0_21" variable_2="k0_21"/>
<!--vars k0_23 and k0_23;-->
<map_variables variable_1="k0_32" variable_2="k0_32"/>
<map_variables variable_1="k0_25" variable_2="k0_25"/>
<map_variables variable_1="k0_52" variable_2="k0_52"/>
<map_variables variable_1="k0_34" variable_2="k0_34"/>
<map_variables variable_1="k0_43" variable_2="k0_43"/>
<map_variables variable_1="k0_45" variable_2="k0_45"/>
<map_variables variable_1="k0_54" variable_2="k0_54"/>
<map_variables variable_1="k0_56" variable_2="k0_56"/>
<map_variables variable_1="k0_65" variable_2="k0_65"/>
<map_variables variable_1="k0_61" variable_2="k0_61"/>
<map_variables variable_1="k0_16" variable_2="k0_16"/>
<map_variables variable_1="delta" variable_2="delta"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<!--vars alpha_pp and alpha_pp;-->
<!--vars C_T and C_T;-->
<map_variables variable_1="n" variable_2="n"/>
<map_variables variable_1="z_c" variable_2="z_c"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="mu" variable_2="mu"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="Cell_concentration"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
<map_variables variable_1="J_CFTR" variable_2="J_CFTR"/>
<map_variables variable_1="J_CFTR_HCO3" variable_2="J_CFTR_HCO3"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="Apical_concentrations"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_CFTR" variable_2="I_CFTR"/>
<!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
</connection>
<connection>
<map_components component_1="ENaC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_ENaC" variable_2="J_ENaC"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_ENaC" variable_2="I_ENaC"/>
</connection>
<connection>
<map_components component_1="IK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_IK" variable_2="J_IK"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Basolateral_concentrations"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_IK" variable_2="I_IK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_BK" variable_2="J_BK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_concentrations"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_BK" variable_2="I_BK"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="GLUT2"/>
<map_variables variable_1="glucose_s" variable_2="glucose_s"/>
<!--vars J_GLUT and J_GLUT;-->
</connection>
<!-- def map between interstitium_concentration and GLUT2 for
vars glucose_s and glucose_s;-->
<!-- //vars J_GLUT and J_GLUT;
enddef;-->
<!-- def map between interstitium_concentration and parameters for
vars time and time;
vars v_cell and v_cell;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="GLUT2"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<!-- def map between GLUT2 and Basol_voltage for
vars I_GLUT and I_GLUT;
enddef;-->
<!-- def map between Cell_concentration and Na_Cl for
vars Na_i and Na_i;
vars Cl_i and Cl_i;
vars J_NaCl and J_NaCl;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="NHE3"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="H_int" variable_2="H_int"/>
<map_variables variable_1="H_Ext" variable_2="H_Ext"/>
<!-- vars HCO3_i and HCO3_i;
vars HCO3_m and HCO3_m;-->
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NHE3_Na" variable_2="J_NHE3_Na"/>
<map_variables variable_1="J_NHE3_H" variable_2="J_NHE3_H"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NaK_channel"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NaK_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
<map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="NaK_channel"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NaK_channel"/>
<map_variables variable_1="I_Na_NaK" variable_2="I_Na_NaK"/>
<map_variables variable_1="I_K_NaK" variable_2="I_K_NaK"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NBC"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
<map_variables variable_1="J_NBC_HCO3" variable_2="J_NBC_HCO3"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NBC"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NBC"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_NBC_HCO3" variable_2="I_NBC_HCO3"/>
<map_variables variable_1="I_NBC_Na" variable_2="I_NBC_Na"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="NBC"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
</connection>
<!-- def map between Cell_concentration and NBC for
vars Na_i and Na_i;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="Apical_Diffusion"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="J_D_Na" variable_2="J_D_Na"/>
<map_variables variable_1="J_D_K" variable_2="J_D_K"/>
<map_variables variable_1="J_D_Cl" variable_2="J_D_Cl"/>
</connection>
<!-- def map between Apical_concentrations and Cell_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Apical_concentrations" component_2="Apical_Diffusion"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Apical_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="P_Na" variable_2="P_Na"/>
<map_variables variable_1="P_k" variable_2="P_k"/>
<map_variables variable_1="P_Cl" variable_2="P_Cl"/>
</connection>
<!-- def map between parameters and Cell_concentration for
vars A_Apical and A_Apical;
enddef;-->
<!-- def map between parameters and Apical_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Apical_Diffusion" component_2="Apical_voltage"/>
<map_variables variable_1="I_D_Na" variable_2="I_D_Na"/>
<map_variables variable_1="I_D_K" variable_2="I_D_K"/>
<map_variables variable_1="I_D_Cl" variable_2="I_D_Cl"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Apical_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="J_D_Bl_Na" variable_2="J_D_Bl_Na"/>
<map_variables variable_1="J_D_Bl_K" variable_2="J_D_Bl_K"/>
<map_variables variable_1="J_D_Bl_Cl" variable_2="J_D_Bl_Cl"/>
</connection>
<!-- def map between Basolateral_concentrations and Cell_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basolateral_concentrations" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Basol" variable_2="A_Basol"/>
<map_variables variable_1="P_Bl_Na" variable_2="P_Bl_Na"/>
<map_variables variable_1="P_Bl_k" variable_2="P_Bl_k"/>
<map_variables variable_1="P_Bl_Cl" variable_2="P_Bl_Cl"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Cell_concentration"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NHE3"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<!-- vars time and time;
vars R and R;
vars T and T;-->
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="NHE3"/>
<!--vars v_sc and v_sc;-->
<map_variables variable_1="I_NHE3_H" variable_2="I_NHE3_H"/>
<map_variables variable_1="I_NHE3_Na" variable_2="I_NHE3_Na"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="NHE3"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<!-- def map between Cell_concentration and NBC for
vars Na_i and Na_i;
enddef;-->
<!-- def map between AE1 and NBC for
vars HCO3_s and HCO3_s;
vars HCO3_i and HCO3_i;
enddef; -->
<connection>
<map_components component_1="parameters" component_2="AE1"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<!--vars A_Apical and A_Apical;-->
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="AE1"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
<map_variables variable_1="J_AE1_Cl" variable_2="J_AE1_Cl"/>
<map_variables variable_1="J_AE1_HCO3" variable_2="J_AE1_HCO3"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="AE1"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<!-- def map between Apical_voltage and AE1 for
//vars v_sc and v_sc;-->
<!-- vars I_AE1_Cl and I_AE1_Cl;
vars I_AE1_HCO3 and I_AE1_HCO3;
enddef;-->
<!-- def map between Basolateral_concentrations and Basol_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basolateral_Diffusion" component_2="Basol_voltage"/>
<map_variables variable_1="I_D_Bl_Na" variable_2="I_D_Bl_Na"/>
<map_variables variable_1="I_D_Bl_K" variable_2="I_D_Bl_K"/>
<map_variables variable_1="I_D_Bl_Cl" variable_2="I_D_Bl_Cl"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Basol_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<!--vars F and F;-->
</connection>
<!-- def map between Apical_concentrations and parameters for
vars time and time;-->
<!-- //vars F and F;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="P_P_Na" variable_2="P_P_Na"/>
<map_variables variable_1="P_P_k" variable_2="P_P_k"/>
<map_variables variable_1="P_P_Cl" variable_2="P_P_Cl"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
</connection>
<!-- def map between Apical_concentrations and Paracellular_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Paracellular_Diffusion" component_2="Paracellular_voltage"/>
<!-- vars I_D_P_Na and I_D_P_Na;
vars I_D_P_K and I_D_P_K;
vars I_D_P_Cl and I_D_P_Cl;-->
<map_variables variable_1="v_ms" variable_2="v_ms"/>
</connection>
<!-- def map between parameters and Paracellular_voltage for
vars capacitance and capacitance;
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<!--vars F and F;-->
</connection>
<!-- def map between interstitium_concentration and Blood_concentration for
//vars J_G_Blood and J_G_Blood;-->
<!-- vars glucose_s and glucose_s;
enddef;-->
<!-- def map between parameters and Blood_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_A_GLUT" variable_2="J_A_GLUT"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
<!--vars J_GLUT and J_GLUT;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="A_GLUT2"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<!-- def map between Cell_concentration and Hydration for
vars H_int and H_int;
vars HCO3_i and HCO3_i;
vars J_buf_C and J_buf_C;
vars H_Ext and H_Ext;
vars HCO3_m and HCO3_m;
vars J_buf_A and J_buf_A;
enddef;-->
<!-- def map between parameters and Hydration for
vars time and time;
vars v_cell and v_cell;
enddef;-->
<!-- def map between Paracellular_Diffusion and Hydration for
vars CO2_i and CO2_i;
vars CO2_m and CO2_m;
vars J_D_P_CO2 and J_D_P_CO2;
enddef;-->
</model>
